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Abstract 

A Brownian motion model is proposed to study parametric correlations in 
the transmission eigenvalues of open ballistic cavities. We find interesting 
universal properties when the eigenvalues are rescaled at the hard edge of the 
spectrum. We derive a formula for the power spectrum of the fluctuations 
of transport observables as a response to an external adiabatic perturbation. 
Our formula correctly recovers the Lorentzian-squared behaviour obtained by 
semiclassical approaches for the correlation function of conductance fluctua- 
tions. 
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I. INTRODUCTION 



The theory of random matrices has been anphed to many branches of physics, such 
as complex nuclei0, disordered metalhc grains&a, the general theory of quantum chaotic 
system^l, disordered conductordl^. random surfaces^, QCdIHI and more recently strongly 
correlated many-particle systemsli3'll3. One of the most striking and robust predictions of 
random matrix theory (RMT), which has far reaching implications in all these different fields 
is the celebrated Wigner-Dyson statistics0'El, which states that the probability distribution 
of level spacings (usually regarded as energy eigenvalues of some random Hamiltonian) is 
a universal function if the eigenvalues are measured in units of the mean level spacing. In 
the case of quantum chaotic systems, it is believed that the success of RMT depends on 
the chaotic dynamics of the equivalent classical system. This dependence seems to be so 
strong that many results derived from random matrix ensembles, such as the Wigner-Dyson 
statistics, have been used as signatures of quantum chaosi. 

An interesting extensionliBEj of RMT is to consider variations in the energy levels of 
the physical system resulting from external adiabatic perturbations. These perturbations 
could be, for instance, a change in the Aharonov-Bohm flux through a ring, or a slight 
alteration in the position of some impurity in a disordered metallic grain, or a variation 
in the strength of an external applied magnetic field or even simply a tiny modification in 
the geometry of the sample. The theory, subsequently developed by Simons, Altshuler, Lee 
and othersij to account for this effect, has proven to be a significant extension of RMT, 
which we shall call parametric random matrix theory (PRMT). Remarkably, it has been 
demonstrated that PRMT, while incorporating all features of RMT, predicts new universal 
statistics after appropriate rescaling of the physical parameters, and thus provides an even 
more powerful characterization of quantum chaotic dynamics. PRMT contains only two 
system dependent parameters: the mean level spacing. A, and the mean square gradient of 
levels, Co, which is defined as the ensemble average 

where Ei = Ei/A denote renormalized energy levels and U is the parameter that controls the 
strength of the external perturbation. A striking prediction of PRMT is that the n-point 
correlation function for density of states fluctuations becomes a universal expression if the 
energy levels are measured in units of A and the perturbation parameter, U, is renormalized 
by Cq. This universality has been derived analytically for disordered metallic grains with 
orthogonal, unitary and symnlectic symmetriesH^^ and has been verified numerically in a 
number of chaotic systemglHlllJH3EirE3. 

We remark that the universality obtained in the framework of PRMT, and similarly that 
of Wigner-Dyson statistics, applies only for levels away from the edge of the support of the 
spectrum. This statement can, in principle, be formally justified by means of renormalization 
group arguments. A simple renormalization group procedure has recently been devised by 
Brezin and Zinn-JustinS. These authors have shown that there are two kinds of fixed points 
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in the renormalization group equations: a stable gaussian one governing the behavior of the 
system at the bulk of the spectrum and an unstable one governing a small crossover region 
around the endpoint of the support where the average density goes to zero as a power 
law. The attractive gaussian fixed point supports the general validity of the Wigner-Dyson 
statistics everywhere at the bulk of the spectrum. It is therefore natural to expect that 
the results of PRMT are of similar general validity, although an explicit proof is not yet 
available. 

At the edge of the spectrum, however, Wigner-Dyson statistics breaks down and an 
entirely new regime of universal statistics emerges. There are two kinds of edges in RMT: 
a hard edge and a soft edge&0, where the eigenvalue support is bounded and unbounded 
respectively. The Gaussian ensembles have two soft edges, while the Laguerre ensembles 
exhibit a hard and a soft edge. There are also ensembles with two hard edges, like for 
instance the Jacobi ensembles. 

The hard edge of the Laguerre ensemble is important in the description of many physical 
systems such as disordered metallic conductors, ballistic cavities and In a recent 

paper Slevin and Nagao0 have introduced a group of matrices, which they called f2-matrices, 
that gives rise to a very powerful way of studying the Laguerre ensemble. Although they 
have used the mathematical structure of this group to propose a model to describe disordered 
metals, we believe that the hard edge of the ensemble generated from the group of f2-matrices 
is actually more appropriate to describe ballistic cavities. 

Ballistic cavities are of considerable current interest mainly because of important recent 
breakthroughs in nanolithograpv. which has enabled its construction in novel high-mobility 
semiconductor heterostructuresE3. These systems have elastic and inelastic mean free paths 
exceeding the device dimensions at sufficiently low temperatures. As a consequence, trans- 
port in these structures is dominated by scattering at the boundaries of the sample rather 
than the more usual impurity scattering of mesoscopic metals. Two striking features of the 
physics of these systems are weak-localizationl^l and universal conductance fiuctuations@@. 
Such typically mesoscopic phenomena appear in ballistic cavities due to the chaotic nature 
of the boundary scattering potential. 

In this work we shall mostly be concerned with a set-up that consists of an interac- 
tion region of finite volume (the resonant cavity in a microwave experiment or the ballistic 
microstructure in mesoscopic physics) connected to two reservoirs by free propagation re- 
gions (wave guides for microwaves or perfectly conducting leads for electron waves) in which 
asymptotic scattering channels can be defined. The interaction region is assumed to "trap" 
the incoming waves by irregular boundary scattering thereby driving the system to a regime 
where the ray optic limit (or classical dynamics) is dominated by classical chaos. We shall 
call such a set-up an open ballistic cavity. In particular, we study the problem of paramet- 
ric correlations in open ballistic cavities, that is, the response of the random spectra to an 
external adiabatic perturbation. Our main objective is twofold: first, we want to extend 
the PRMT to describe open quantum chaotic systems, such as ballistic cavities coupled to 
external reservoirs; and second we want to understand how the theory changes when the 
region of physical interest is at the hard edgeof the support of the spectrum. The physics 



of the soft edge has been discussed elsewhereca. 

Parametric correlations in open mesoscopic systems has recently been the subject of 
many works&0. A concise account of some of our results has been presented previouslyi^. 



3 



In this work we provide much more details of the calculations, a different interpretation of 
the hydrodynamic limit and new results. 

In section II we discuss the S-matrix approach to ballistic systems. In section III, we 
review some of the properties of the group of f2-matrices of Slevin and Nagao. In addition, 
we derive a diffusion equation generated by a random walk on the fi-matrix manifold. We 
propose a modified version of this equation as a Brownian motion model for parametric 
correlations in open ballistic cavities. In section IV, the Brownian motion model is used to 
derive exact non-perturbative expressions for the parameter-dependent two-point correlation 
function at the hard edge of the spectrum. We demonstrate that after appropriate rescaling 
this function becomes universal. In section V we show that in the hydrodynamic limit our 
theory predicts a Lorentzian-squared behaviour for the correlator of conductance fluctuations 
of ballistic cavities in agreement with semiclassical calculations. Finally, we derive a formula 
for the power spectrum of an arbitrary linear statistic on the transmission eigenvalues. A 
summary and conclusions are presented in section VI. 



II. S-MATRIX AND BALLISTIC SYSTEMS 

It is now well established that the most convenient description of quantum transport in 
open ballistic cavities is given by the S-matrix ^30. By deflnition the ^-matrix relates the 
incoming flux amplitudes /; and to the outgoing ones Oi and Or through the formula 

s(i-)-("'), (2) 



where the subscripts / and r denote the left and the right sides of the sample respectively. 
Current conservation implies that S is unitary: = S'' . A general and rather simple 
explicit parametrization for S isi^ 



where r and r' are N x N reflection matrices; t and t' are N x N transmission matrices; 
u^^\i = 1, . . . , 4) are N x N unitary matrices and r denotes an x diagonal matrix with 
real eigenvalues < Tq, < l(a = 1, 2, . . . , A^), which are called transmission eigenvalues. The 
total transmission coefficient can be written as 

N 

T = trttt = ^ r^. (4) 

a=l 

Therefore, the Landauer-Biittiker conductance is simply 

N 

G = Go E (5) 

a=l 

The remarkable simplicity of this formula, in comparison with the Kubo-Greenwood 
expressionii for the conductance obtained directly from linear response theory, highlights 
the advantage of the Landauer-Biittiker formalism. This, however, cannot be the end of the 
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story because in principle the transmission eigenvalues, Tq,, depend in a very complicated way 
on the underlying quantum dynamics induced by the scattering mechanism. Fortunately, in 
some cases of interest such as ballistic cavities and quasi-one-dimensional disordered con- 
ductors the complexity of the dynamics turns out to be such that most details of microscopic 
origin are irrelevant and the joint probability distribution of the t^s can be determined from 
symmetry arguments which are well known in the context of random matrix theory!!. 

In the case of ballistic cavities, it has been shownS that the S'-matrix fluctuates as a 
function of the incident momentum due to multiple overlapping resonances in the cavity. 
There is considerable numerical evidenceil^ to support the assumption that through this 
fluctuations the S'-matrix covers its manifold with uniform probability. With this ergodic 
hypothesis in mind, we can study the fluctuations in the S'-matrix by means of an ensem- 
ble of random matrices with a uniform distribution. So, dPj3{S), which is the probability 
to find S in the volume element dS, is proportional to the Haar measureS of the group 
of S-matrices. In terms of the parametrization (^) one finds, after integration over the 
eigenvectors distributionsii'0 

dP^ir) = dfipir) = CnMt) n ri^-'^/' 11 dra, (6) 

c a 

where Cn^js is a normalization constant and (3 is a symmetry index, whose value is /5 = 1 for 
systems with time-reversal symmetry (T-symmetry) in the absence of spin-orbit scattering, 
(5 = 2 for systems without T-symmetry and (3 = A for systems with T-symmetry in the 
presence of spin-orbit scattering. The factor Jpir) = Y[a<b ka — is ultimately responsable 
for transmission eigenvalue repulsion in the ensemble. 

At this stage it is convenient to make the following change of variables, = 1/(1 + uf), 
with < Ui < oo. From Eq. @ we find the following joint probability density for the 
variables {ui} 

P{iy) = Z-^exp{-f3n), (7) 

where 

^ = -^E<5(^-^.) + En^0, (8) 



Q{u,u') = ln\u^ -u% (9) 

and 
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V{iy) =[N + j ln(l + - ^ In ly. (10) 

Note that P{iy) has the form of a Gibbs distribution and Ti plays the role of a Hamiltonian 
of classical particles with logarithmic pairwise repulsion, Q(z/, z/'), and a confining potential, 
V^(z/). The dimensionless Landauer-Biittiker conductance in this new set of variables reads 

N 

9 = T.9{^i)^ (11) 

i=l 
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where (?(z/) = (l+z/^)~^. In RMT observables of this form are called linear statistics^ because 
products of different eigenvalues do not appear in their defining expressions. Consequently, 
their statistical properties are the simplest. In particular, the average and variance of Eq. 
(|ll|) are simply 

(g) = / g{u)p{u)du, (12) 

JO 



vaiig) = g{u)g{u')S{u,u')dudiy', (13) 



in which p(z/) and Slu, u') are the average level density and two-point correlation function 
respectively. There is a very powerful way of calculating ^(z/, z/') which has been developed 
by Beenakkeril. It is based on the following identity 

which can be easily verified. For ^ 1, one can show that 

V{u) = / p{u')Q{iy,u')diy'. (15) 
Jo 

Using ([1^) the functional derivatives in (|1^ can be performed and we find 

V + u'' 



u — u 



+ 0{1/N). (16) 



A remarkable consequence of this technique is the conclusion that the leading expression 
for S{v,v'), in an expansion in inverse powers of A^, is independent of the potential V{v). 
Therefore, the variance of a linear statistic, such as the conductance (see Eq. (p^)), is the 
same for all random matrix ensembles with an edge at the origin of the spectrum (since 
the i/'s are all non-negative) and with an eigenvalue repulsion potential equal to Q{i>,iy'). 
Consequently, if we are only interested in studying the dominant contribution to fiuctuations 
in mesoscopic observables, we are free to choose the functional dependence of Vlu) that we 
find convenient. It turns out that the simplest choice for V{i') can be obtained from (p!0D 
by means of the rescaling, u — > u/N, and by taking the large A^ limit keeping u fixed. We 
then get 

V{u) = ^u'-^\nu. (17) 

We remark that Eqs. (^, (|) and (^ with V(z/) given by ([17D constitute the Laguerre 
ensemble of random matrices^. This ensemble has a number of interesting universal prop- 
erties which will be discussed in section III, where it is used to build a model of parametric 
correlations in ballistic cavities. 

We conclude this section by giving the leading terms in the expansions of (g) and vai{g) 
in inverse powers of A^ 
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{9) = \n+^ + 0{1/NI (18) 

var(^?) = ^ + 0(l/iV). (19) 

The leading term in Eq. ( [T^ ) is a direct consequence of the ergodic hypothesis, since it 
imphes that the average transmission and reflection probabilities are equal. The second term 
in Eq. (18) corresponds to the weak- localization correction, which as is due to constructive 



interference of time-reversed backscattering trajectories. Eq. (^) is an illustration of the 
remarkable phenomenon of universal conductance fluctuations, which is a clear signature of 
the non-self-averaging nature of observables in quantum transport theory. 



III. THE GROUP OF ^-MATRICES 

In this section we introduce the group of f2-matrices as an abstract mathematical struc- 
ture, which, as we shall see, is very convenient for studying the Laguerre ensemble. We 
remark that our approach differs considerably in philosophy from that used by Slevin and 
Nagao0. They have motivated the group of r2-matrices by relating it to the group of trans- 
fer matrices of a disordered conductor. We believe that to make such a connection from 
the outset is unnecessary and might lead to misinterpretations. The group of f2-matrices 
is interesting in its own right as a powerful way of describing statistical properties of the 
Laguerre ensemble, which we regard as a tractable mathematical model for studying local 
eigenvalue correlations in systems with a hard edge in the spectrum. Furthermore, as we 
discussed in section II, the formula for the variance of linear statistics derived from the 
Laguerre ensemble applies quite generally to any maximum-entropy ensemble with a hard 
edge, most notably, the ensemble describing ballistic cavities. 

By definition, a general complex matrix VL satisfies 

S^l^S^ = -n. (20) 
If the system has T-symmetry we have the additional relation 

VL* = S^l^S^. (21) 



The matrices and S^. are given by 



and 



E.= (« I), (23) 

where 1 is the N x N unit matrix for orthogonal and unitary systems, and it is the N x N 
unit quaternion matrix for symplectic ensembles. 
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One can show that a general matrix Q satisfying ( pOD and (pi]) can be written in the 
form 











)(° 




[ 





"=(o ojio 

where u and f are N x N unitary matrices for unitary and orthogonal systems and they 
are N x N quaternion unitary matrices for symplectic systems. The matrix u is a. diagonal 
matrix with real eigenvalues. One can easily check that the matrices Q form a group under 
addition, which we call Q^{fl, +). 

Let fit denote a point on the f2-matrix manifold that is generated from a random walk 
defined by 

nt = T.^s{st), (25) 

s=l 

where St is the step of the walk and Qs{St) is a random matrix belonging to Q^{Q, +). The 
simplest choice for Qs{St) is 

^],(5t) = 2^v^(_°^ ly (26) 

where y is a complex random matrix with vanishing average and with second moment given 
by 

1 

{yabVcd) = ^(^jy ^ i^'^cSbd + SadSbc), = 1, (27) 



{yabVcd) = -^^aAd, 



(3 = 2, 



(28) 



{VabaVcda') 



a,a' 



(29) 



in which 7 is a constant. Note that for systems with symplectic symmetry {(3 = 4) we have 
used the fact that y is a quaternion matrix and thus any of its matrix elements can be 
written as yab = J2a=oyabaGa, where are the quaternion units. 

We can rewrite (p5|) as Qt+st = flt + fls{St), which in view of ( p4D implies, to order 0{6t), 
the stochastic process 



Sui = iV5t{Au - A\)lv, + 



2bt 



7(/3iV + 2-/3) 



1 



Vi 



(30) 



where Aij = {u^yvh')ij and z/j = {y)ii. If we take the continuum limit of (pOf ) we find the 
following Fokker-Planck equation for the probability distribution Wlu, r) 
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dr 

where r = 2t/{PN + 2-13) and 



7?^ = E(-7^-A + £,)H^, (31) 




A = - + 2/3 5] (32) 
It is interesting to note that the transformation W = JP imphes that P satisfies 

where J = Ylii'iYliKj Wf ~ ^]\^ and is the Laplace-Beltrami operator of the fi-matrix 
manifold. 

In the analysis of Slevin and Nagao0, the distribution of ^-matrices was determined by 
a maximum-entropy criterion to be the Laguerre ensemble 

W^v) = C^,, n - ^if n ^t^" ^M~cu^). (34) 

i<j i=l 

where Cat^/j and c are constants, a = and P is the usual symmetry index. This result 
suggests, in analogy with transport theory, that our Brownian motion in the Q manifold 
must satisfy the requirement that when r — >■ oo, W{v, r) — > W^{v), that is the distribution 
function evolves to a configuration that maximizes the entropy. This can be achieved simply 
by adding to Di a term corresponding to a confining force so that (^) becomes 

dW d / ldW\ 

in which 

m = -(« + ^) E In + ^ E - E In - v] I . (36) 

where (3 = 1 and a = 0. Eq. (|35| ) is the main result of this section. It can be regarded 
as representing a Brownian motion of classical particles at positions {^'1(1")} moving in a 
viscous fiuid with friction coefficient 7"^ at temperature (3~^. The parameter r is a fictitious 
time which we shall relate to the perturbation parameter, so that (|35| ) becomes a model to 



describe parametric correlations in the Laguerre ensemble. In the framework of PRMT it 
has been explicitly demonstrated that all the universal functions can be obtained from a 
Brownian motion model if one assumes that r = (5M^)il^, where 5u is a parameter charac- 
terizing the strength of the external perturbation. Motivated by this result, we assume the 
same relation to apply to our model. We stress that the consistency of this assumption has 
recently been confirmed by explicit calculations of a similar model using supersymmetryS. 
The apparently artificial parameter a has been introduced to allow the classification and 
understanding of the universal expressions, which we derive in section IV, in the more gen- 
eral context of Laguerre ensembles. Note that Eq. (|3^) with a > — 1 characterizes the most 
general Laguerre ensemble of random matrices. 
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IV. UNIVERSAL CORRELATIONS AT THE HARD EDGE 



A crucial difference between open ballistic systems and metallic conductors is that, in 
the former, the joint probability distribution of the transmission eigenvalues can be directly 



obtained from a maximum-entropy principle, whilst as shown in Refs. ^ and ^ 

this is not the case for disordered conductors. A direct consequence of this fact is that the 
fluctuations in transport observables of ballistic systems are completely characterized by the 
universal logarithmic repulsion of conventional random matrix ensembles, like the Laguerre 
one. An outstanding common feature of the ensembles of random matrices appropriate to 
describe open ballistic systems and disordered conductors, though, is that the spectrum of 
transmission eigenvalues has a hard edge, since the eigenvalues are all non-negative. This is, 
as we discussed in the Introduction, also a feature of the Laguerre ensemble. In this section 
we shall demonstrate that the presence of this hard edge implies the existence of universal 
parametric correlations of a new kind. 

A. Exact Mapping onto a Schrodinger Equation 

Our starting point is Eq. (p5D-(^), which we rewrite here in the form 




dW 

7^ = CfpW, (37) 
where Lpp is the Fokker-Planck operator, given by 

(38) 

i UAi \ UAi J 

where = z/f and 

D W = 4 5] ^\ - 2c\i + 2a + A/ 13 (39) 

while 

Df> = fs,. (40) 
At this stage, it is useful to perform the following transformation^ 

P(A,r)=exp(-i$(A))v^(A,r), (41) 
which implies that \I'(A,r) satisfies the evolution equation 

^ = C^, (42) 

where £ = £Ms given by 
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L = e*/2£i.pe-*/^ (43) 
From (|3|), (|D-(|0|) and (||) we obtain C = B - H, where 

5 = Nc{2 + (3{N + a - l))/(27), uo = Ac/ {(3'^) and = cAj. The operator H can be 
regarded as a one dimensional Hamiltonian of fictitious hard-core particles, for which it is 
convenient to assign fermionic statistics. We remark that it is a general property of the 
mapping (^) that equilibrium fluctuations of the classical Brownian particles correspond 
to quantum mechanical ground-state fluctuations of these fictitious fermions. Note that for 
f3 = 2 the two-body interaction term vanishes and Ti becomes the Hamiltonian of a system 
of free fermions. In this case we can rewrite Ti. as 



oo 



n = J2 ^p4^p^ (45) 

p=0 

where Ep = (2c/7)(j9 + (a + l)/2) and cj^ (c^) creates (annihilates) a fermion with quantum 
number p. It is convenient at this stage to introduce the field operators and 
where ip{X) = Y.p4>p{^)Cp and = Y.p'PpWcl, in which 

with Lp{x) denoting the associate Laguerre polynomial. 

With this notion, we turn to the calculation of some quantities of physical interest, 
namely the average level density and the two-point correlation function for level-density 
fluctuations. In second quantized language the local density operator is defined as 

hiX,T)=^P\X,T)^lj{X,T) (47) 

where 

ij\X,T) = e^^ij\X)e-^^ (48) 

and 

^(A,r) = e^"V(A)e-^". (49) 
Let |0) denote the A^-fermion ground-state of (^Sf ). Then the average level density is simply 



p(A,r) = (0|n(A,r)|0)= E0,(A)2 = p(A), (50) 

p=0 

which is independent of r, as expected, since the Brownian particles are in a stationary 
state. 
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The two-point function for level-density fluctuations is defined as 

^(A, A', r) ^ (0|n(A, r)n(A', 0)|0) - p(A)p(A). (51) 
Using (|47| ) and Wick's theorem we find 



^(A, A', r) = Go(A, A', r) E.to' 0,(A)</),(A')e^'- - 

EplV EJTo' 0p(A)0,(A')0,(A)0,(A')e(^^-^)^ (52) 



where 



w ,_ cexp(-f(A + A0coth(.;r/2)) ^ f ^ 

^ ' - 2sinh(a;r/2) Uinh(^^/2) j ^ ^ 



and /a (a;) is a Bessel function. Equations ( E^ ) and ( pB] ) can be considered as a generalization 
to finite r of the standard result of RMTS, which we recover by taking the limit r — > 0, 
giving 



^(A, A', 0) = 6{\ - A')p(A) - {k{\, A')) , (54) 

where 

N-l _ 

ki\,\') = Y.4>pWM>^'). (55) 

p=0 

In the next section we show how to rescale the average level density and the two-point 
correlation function at the hard edge of the spectrum. 

B. Rescaling Physical Quantities at the Hard Edge 

In order to obtain the new universal behaviour anticipated in the Introduction we need 
to rescale all the quantities of physical interest at the hard edge of the spectrum, that is at 
A = 0. This is done more easily if we reintroduce the original variables z/j = a/A^. So, the 
average level density becomes 

7V-1 

piu) = y: <pm\ (56) 

p=0 

where (ppi^v) = v^2z/0(z/^). For large we can use the asymptotic^ 

(t>N{y) ^ {2cuf'^J^ (2(iVc)i/V) , (57) 

where Jq,(x) is a Bessel function. Inserting (|57|) into (|56|) and taking N ^ oo and c — 
such that 2y/Nc — >• pon we find 

p(z/) = vrVo / dssuj^npous), (58) 
Jo 



12 



which is exactly the result obtained in Refs.ETEj for the average level density rescaled at the 
hard edge of the spectrum. Note that the constant po is, as can be seen from the relation 
~ pq, valid for large z/, just the bulk average level density. 
A similar calculation yields for the two-point function 

S{iy, iy',6u) = tt^PqJq ds ds'ss'i'iy'Ja{iTpoi's)Ja{iTpoi''s)Ja{TipQiys') 

Ut^Po^'s') exp [^6u\s' - s")) . (59) 



This is the central result of this section. Note that, for Su = 0, Eq. ( |5D| ) reproduces 
known resultsiHil for the two-point correlator of ensembles with a hard edge. Therefore, 
we believe that Eq. (|59D represents an extension of these RMT results to account for the 
dispersion of the eigenvalues as a function of the perturbation parameter. Bearing in mind 
physical applications, consider now the particular case of the ensemble describing open 
ballistic sytems, which can be mapped (see section 2.2) onto a Laguerre ensemble with 
a = 0. Following Ref. we make the rescalings z> = p^u, v' = pou', 5u^ = du^pl/'j and 
S^O, z>', 6u) = Si^u, u', 5u) I Po- We can see that S{y, z>', bu) is a function that is independent of 
any physical parameter, such as the Fermi velocity or size of the sample, and therefore can be 
regarded as a universal characterization of quantum chaotic scattering in such systems. For 
the more general case of Laguerre ensembles, with a > —1, we can see that the parameter 
a labels new universality classes characteristic of systems with a hard edge in the spectrum. 
Finally, we remark that Eq. ( |59D has recently been derived, for a = 0, by Andreev, Simons 
and Taniguchiii using the supersymmetry technique. In addition, they have demonstrated 
that the universal behaviour at the hard edge extends to another kind of two-point function. 
Another quantity of some interest is the distribution of level velocities defined as 

K{v) = {6{v-du/du)). (60) 

One can show thatS 

K(v) = lim , . Siu,!/ + v6u,6u). (61) 



Now using ( p2D with r = 5u we find 



S{v, V + vbu, bu) ^ 4^ ^'^ exp f-^) (62) 



bu \2tx J 

and therefore the distribution of level velocities is gaussian 

^ \l/2 / ryy'i^ 



Remarkably, this is exactly the same result that one finds at the bulk of the spectrum. 
We stress that this result, although simple, is by no means trivial, since the hard edge 
rescaling procedure introduces considerable changes in most physical quantities, such as the 
average level density and two-point correlation function. One interesting application of (|63[) 
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is associated with the interpretation of the rescahngs discussed right after Eq. (|59D. The 
quantity Cq = Po/7 can now, in view of (^), be written as 

Co = pI{v'), (64) 

so that the rescahngs = pqu and Su"^ = Sv^Cq mean that we get universal functions if 
we measure the levels in units of the average bulk level spacing p^^ and if the perturbation 
parameter is rescaled by the mean square gradients of the levels \/Cq. 



V. THE POWER SPECTRUM FORMULA 



A complete characterization of the fluctuations in transport observables in mesoscopic 
systems consists of a description of both their magnitude and power spectrum. One can get 
this information from the correlator 



F{5u) = {SA{Su)SA{0)), 



(65) 



where A is an arbitrary transport observable and 6u parametrizes the external perturbation 
on the system which drives the fluctuations. The power spectrum C{uj) is defined as 



C{UJ) 



dte"^'F{t). 



(66) 



Microscopic diagrammatic calculations@@ for disordered metals suggest that F{5u) decays 
as a power law as a function of 5u/£c-, where £c is a correlation parameter which sets the 
typical scale of the spacings between the peaks and valleys in the values of A, for a given 
realization, as a function of 5u. This slow power law decay of F{6u) indicates that the 
system sustains some sort of long-range memory in the fluctuations of the observable. This 
kind of behavior must be contrasted with the exponential decay of Poisson processes and 
can be regarded as a manifestation of spectral rigidity, a well known phenomenon in the 
theory of random matrices. 

It is well known@'^ that universal fluctuations in transport observables of mesoscopic 
disordered conductors can be obtained by simply taking the hydrodynamic limit of the 
DMPK equation. We expect a similar procedure to apply to the power spectrum formula of 
open ballistic structures. In section II, we have demonstrated that the Laguerre ensemble 
can be used to describe fluctuations in transport observables of ballistic cavities. With this 
in mind we shall now discuss the hydrodynamic limit of Eq. (|35D . 

Multiplying (E^j by J2i ^{^ ~ ^i), integrating over all z/j's and using the definition 



(67) 



where (. . .)r denotes an average over the distribution W{h',T), we get in the large N limit 
the following evolution equation for the average level density 



7- 



Or 



d_ 
du 



p{u,T 



d_ 
du 



dpp{p, t) In 



2 2 




v — p 


-¥)] 



(68) 
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One can see from (|68|) that the equihbrium density satisfies 



Peg(/i) In 



2 2 
U — fi 



C 2 

dfi = -u + const. 



(69) 



The constant on the right hand side of (^) is determined by the normahzation condition 

Pe,{u)du = N (70) 



One can easily verify that 

' 1 for z/ > 2(iV/c)i/2 ^ ^ 

is the solution of ( |BUD satisfying ([7DD. Note that for large N at fixed one finds 

' ^ (72) 



Peg{u) ~ -VA^C = po, 



as expected. In the regime of universal mesoscopic fiuctuations that we are concerned with, 
one can safely consider the average level density, po, to be much larger than the fiuctuating 
part (5p(z/, r). Therefore, it makes sense to try and linearize Eq. (|68D by writing p(z/, r) as 



p{u,t) = Po + Sp{u,T). 



Inserting (^3]) into (|68D yields 

d6p{u, r) 



^2 ^oo 



7- 



Since 



dr 



In 



-po 



dv'^ Jo 



2 2 

u — p 



dpSp{p, t) In 



2 / — cos ku cos /cp, 
JO k 



2 2 

z/ — p 



Eq. (|7^ can be solved by means of Fourier cosine transform. We find 

2 /■°° 

5p(z/,r) 



TC JO 



(73) 



(74) 



(75) 



(76) 



dk5p{k, 0) exp(— 7rA;po''"/7) cos ku 

The two-point function S'(i/, z/', 5m) can be obtained from (|76D through the identity 

S{u,u',6u) = {6p{u,T)6p{u',0))^^ (77) 



where (. . .)eq stands for an average over the equilibrium distribution (|^). 
It is useful to define the following double Fourier cosine transform 



S{k, 5u) 



dv / du' S{u, u' , Su) cos ku cos ku' 



(78) 



so that we get from (fT^) 
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S{k, 5u) = (5p{k, 5u^)5p{k, 0)) = exp{-TTkpo5uy-f)S{k, 0) (79) 

\ I eg 

The Fourier component S{k, 0) can be obtained directly from Eq. (|r^), which together with 
(0) yields S{k, 0) = k/2(3, so that 

k 

S{k,6u) = —exp{--Kkpo6u'^/'y). (80) 

Let A = J2i CL{^i) denote an arbitrary observable, which can be expressed as a linear statistic. 
Then from ( |65[ ) and ([77| ) we get for the correlation function 

roc roc 

F{5u) = j^ du dv'a{v)a{v')S{v,v\5u), (81) 

which by virtue of (^Dj) yields the formula 

2 r'^ 

F(5u) = --^ dkkexp{-TTkpo6u^/-f)d^{k), (82) 



where 



a{k) = I coskva{'u)dv. (83) 
Jo 

Finally, from (P^) we can obtain the power spectrum formula 



Equations ( ^2|) and (|84 ) are the principal result of this section. 

As an application, which also serves as a test of our results, we consider the two-probe 
Landauer-Biittiker dimensionless conductance (Eq. (pi])) 

We find 

~a,{k) = ^e-^ (86) 

thus the correlator of Eq. (^) gives 

F{Su) = ^{l + SuyS^y\ (87) 



where 8c = y27/7rpo is the correlation parameter. For the power spectrum we get 

C{u) = ^e-Hl + ujSc). (88) 

We would like to stress that Eqs. (|^) and (|88D, for /? = 2, are in complete agreement with 
independent calculations based on semiclassical quantization!^. 



16 



VI. SUMMARY AND CONCLUSIONS 



In this work we have studied the effects of an external adiabatic perturbation on the 
transmission eigenvalue correlations of open ballistic cavities. 

In particular, we have derived a Brownian motion model to describe dynamic fluctuations 
in the Laguerre ensemble of random matrices, which we have proposed as a model for 
parametric correlations of transmission eigenvalues in ballistic cavities . This model has 
enabled us to obtain explicit non-perturbative expressions for the two-point function of 
level density fluctuations at the hard edge of the spectrum. We have shown that after 
appropriate rescaling this function becomes system independent and can therefore be used 
as a signature of quantum chaos in open ballistic cavities. In the hydrodynamic limit, 
we have demonstrated that the two-point parametric correlator of mesoscopic conductance 
fluctuations in ballistic cavities is a Lorentzian-squared, in agreement with semiclassical 
calculations. We have also obtained a formula for the power spectrum of the fluctuations of 
an arbitrary linear statistic in such systems. 

Parametric correlations in the S-matrix ensemble has recently been discussed in Refs. |7 
and They obtained exponential decay for the conductance correlator, which disagrees 
with equation (|87|) and consequently with microscopic semiclassical calculations. It is not 
understood at the moment why the S-matrix approach, which was so successful in describing 
non-parametric fluctuations, should have a hydrodynamical limit for parametric correlations 
that disagrees with semiclassical results. We remark that our fi-matrix approach, on the 
other hand, although phenomenological in nature, has the advantage of containing the cor- 
rect semiclassical limit. 

The author would like to thank M. D. Coutinho-Filho for reading the manuscript and 
helpful comments. This research was partially supported by the Brazilian Agency CNPq. 
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